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ABSTRACT 

Single  fringe  interferograms  are  presented  of  4-Diil  Cu  wires 
exploded  at  20  kv  into  argon  at  ambient  pressures  of  I/8,  1/16  and 
1/32  atm.  Features  discernible  include  a  compressive  head  shock 
wave,  arc  plasma,  a  weak  plasma  wave  and  the  expanding  metal  wire. 
On  the  basis  of  certain  plausible  assumptions  it  is  seen  that  the 
arc  jjlasma  has  a  temperature  of  about  2.5  ev;  but  its  leading  edge, 
a  region  not  in  thermal  eq,uilibrium,  has  electron  temperatures 

O 

^  lO""  ev  and  is  the  boundary  of  an  electron  driven  shock  wave. 
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INTRODUCTION 


In  a  recent  application  of  interferometric  techniques  f  the 
formation  of  shock  waves  hy  exploding  wires  at  low  ambient  densities 
of  argon,  c.  l/l6  atm,  is  shown  to  depend  on  the  formation  and 
expansion  of  an  arc  plasma  concentric  with  the  wire.  The  interferograms 
presented  there  are  all  of  the  multi-fringe  type  most  useful  for  the 
existing  schemes  of  quantitative  reduction.  Certain  features  of  the 
flow  are  readily  discerned  in  a  multiple  fringe  picture,  but  others 
remain  hidden  from  view  because  of  the  discontinuous  mapping  of  the* 
disturbance  offered  by  the  light  or  dark  fringes.  By  a  slight  modifi¬ 
cation  of  technique  to  employ  a  single  fringe  field,  a  different 
representation  of  the  flow  can  be  achieved  which  exhibits  certain  finer 
details  not  easily  detectable  in  the  multiple  fringe  pictures.  The 
present  study  is  concerned  with  a  semi-quantitative  analysis  of  single¬ 
fringe  interferograms  from  wire  explosions  in  argon,  teiken  under  condi¬ 
tions  corresponding  to  those  of  the  multiple  fringe  interferograms 
presented  in  oin*  previous  paper. 


*  This  paper  will  hereafter  be  designated  by  w- 
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2.  EXPERIMENTAL  DATA 


2 . 1  Int  erf  er  o  grams . 

The  general  technique  is  described  in  [1].  To  obtain  a  single 
fringe  interferogram  one  need  only  increase  the  fringe  spacing  until 
a  single  light  or  dark  fringe,  preferably  one  of  the  central  fringes 
in  the  pattern,  uniformly  covers  the  entire  field  of  the  interfero¬ 
meter.  Plate  imperfections  may  prevent  perfectly  vmiform  illumination 
of  the  field,  but  are  evident  mainly  at  the  edges  and  need  not  concena 
us  when  the  disturbance  is  reasonably  well  centered.  With  these  condi¬ 
tions  of  adjustment  the  interferometer  shows  directly  the  contours  of 
constant  fringe  shift,  and  represents  the  boundary  of  the  disturbance 
in  a  nearly  continuous  curve.  Single  fringe  interferograms  of  4-mil  Cu 
wires  at  20  kv  exploded  in  argon  at  l/52,  l/l6  and  l/8  atm  are  shown 
in  Figures  1-5. 

2.2  General  Features. 

Inspection  of  these  figures  and  comparison  with  the  earlier 
multifringe  pictures  Cl]  allows  certain  major  phenomena  to  be  readily 
identified. 

The  rounded  tip,  foremost  in  the  streak  interferogram,  is  clearly 
the  boundary  of  the  peripheral  arc  within  which  the  negative  fringe 
shift  achieves  its  largest  magnitude.  The  zero  fringe  contour  can  be 
seen  entering  the  pattern  just  outside  the  arc  boundary.  It  is  continu¬ 
ously  connected  to  the  undisturbed  fringe  ahead  of  the  explosion  and 
separates  the  major  area  of  negative  fringe  shift  near  the  axis  from 
the  area  of  positive  fringe  shift  further  out  and  just  behind  the  head 
shock  wave.  The  head  shock  is  the  nearly  straight,  outer  boundary  of 
the  disturbance. 

Numbering  light  fringes  with  half  integers  and  dark  with  whole, 
it  is  clear  that  the  contour  5  =  +  0.5  lies  just  inside  the  shock  on 
each  side  of  the  central  axis.  Then  comes  the  6=0  contoxrr  which 
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originates  in  the  landisturbed  field  ahead  of  the  explosion,  proceeds 
back  and  slightly  outward  only  to  loop  forward  again  at  smaller  radii 
into  the  compressed  region  Just  outside  the  expanding  wire  material. 

In  the  1/32-atm  case  the  returning  zeroth  fringe  occupies  a  large, 
wedge  shaped  area  closest  to  the  weakly  expanding  wire.  Inside  the 
zeroth  contour  lies  that  for  8  =  -  0.5  which  outlines  the  arc,  streajns 
backward  in  time  through  a  broadening  area  and  then  returns  very  sharply 
near  the  expanding  wire  to  connect  with  itself  eind  form  a  closed  loop  at 
the  tip.  Within  it  one  finds  the  5  =  -  1.0  and  8  =  -  1.5  areas  both  of 
which  lie  close  to  the  axis  in  the  region  of  most  intense  luminosity. 

In  the  l/8-atm  case,  luminosity  partly  obscures  the  inner,  negative 
fringe  contours . 

Near  the  tip  of  the  disturbance,  transverse  paths  of  the  inter¬ 
ferometer  beam  are  small  but  negative  fringe  shifts  are  large.  These 
facts  imply  high  electron  densities.  At  large  radii  and  later  times, 
the  paths  are  longer  and  smaller  electron  densities  will  achieve  the 
same  fringe  shift.  The  contour  for  8  =  -  0.5  may  fairly  be  taken  as 
the  area  depicting  the  onset  and  decay  of  the  electron  cloud.  The 
zeroth  contour  represents  the  area  in  which  refraction  effects  of 
electrons  versus  those  of  compressed  ions  and  atoms  just  cancel. 

Two  compressive  shocks  can  be  seen.  For  orientation  refer  to 
Figure  4.  The  first,  which  we  shall  call  SI,  is  the  head  shock  wave 
already  noted.  The  second,  SII,  is  delineated  by  the  sharp  boundary 
near  the  expanding  wire  material  where  fringe  shift  is  positive  in 
passing  from  8=0  to  6=  -0.5  with  increasing  radius.  Inside  the 
shock  SII,  the  wire  material  is  seen  to  expand,  weakly  in  the  l/32-atm 
case  and  more  strongly  in  the  l/l6-  and  l/8-atm  cases.  In  these  latter, 
fringes  of  rather  definite  structure  can  be  seen  within  the  metal 
boiondary  itself. 

Finally,  we  note  a  weak  disturbance,  denoted  by  WS,  originating 
at  the  tip  where  the  8  =  -  0.5  contour  first  closely  approaches  the 
expanding  wire  (see  Figure  2,  the  l/l6-atm  case  particularly). 
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This  weak  disturbance  appears  to  reflect  from  the  wire  and  rise  steeply 
along  a  straight  path,  causing  an  outward  flare  in  the  6  =  -  0.5  contour 
at  about  t  =  2  |J.sec.  A  little  consideration  shows  that  the  fringe  shift 
across  this  disturbance  is  negative,  indicating  a  detectable  increase  in 
electron  density. 

2.5  Measured  Velocities. 

In  the  i)articular  case  of  the  l/l6-atm  interferogram , velocity 
measirrements  have  been  made  of  each  of  the  visible  disturbances  just 

ft 

cataloged  and  are  presented  in  Table  I,  where  P  stands  for  the 


TABLE  I  -  l/l6  ATM  AEGON 


Disturbance 

P 

SI 

SII 

ws 

Velocity  (ram/psec) 

>  20 

1.5 

1.0 

k.2 

Mach  No.  (a^  =  .52  mm/psec) 

>  6o 

plasm  boundary  near  the  tip.  Mach  numbers  have  been  calculated  relative 
to  room  temperature  argon  for  purposes  of  comparison  with  shock  tube 
data.  The  velocity  of  20  ram/psec  for  the  plasma  boundary  P  is  measured 
where  the  front  intersects  the  expanding  wire. 

2,h-  Estimated  Energies. 

Although  it  is  shown  in  our  previous  study  [1]  that  blast  wave 
theory  is  not  strictly  applicable  to  exploding  wire  shocks,  it  never¬ 
theless  provides  a  useful  way  of  classifying  such  shock  waves.  Accord¬ 
ingly  we  give  here  in  Table  II  a  list  of  apparent  energies  in  joules/cm 
indicated  by  the  sequence  of  wire  explosions  into  argon,  reali^iing  that 
the  values  given  do  not  accurately  represent  the  energy  in  the  flow 
field  behind  the  shock.  They  do  have  the  heuristic  merit  of  ordering 
the  shocks  according  to  apparent  energy.  On  the  basts  of  our  previous 
study  it  seems  likely  that  this  ordering  has  validity  for  p  >  l/k  atm. 


TABLE  II 


p  (atm) 

1 

1/2 

1/4 

1/8 

l/l6 

E  (joules/cm) 

6.2 

2.9 

1 

O.T 

0.5 
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5 .  ANALYSIS 


^.1  Plasma  Boundaries. 

We  have  already  noticed  that  the  5  =  -  0.5  fringe  gives  a  crude 
indication  of  the  extent  of  the  electron  cloud  in  any  of  the  single 
fringe  interf erograms .  This  notion  may  be  made  somewhat  more  precise 
by  observing  that  the  zeroth  fringe  contour  is  the  locus  for  which 
the  total  path  difference,  including  negative  contributions  by  electrons 
and  positive  by  compressed  ions  and  atoms,  is  zero.  Near  the  tip  of  the 
distiirbance  where  the  zeroth  fringe  is  first  seen,  the  entire  light  path 
in  the  disturbance  cannot  be  distinguished  from  a  path  of  eq.ual  length 
in  the  \indisturbed  beam.  Therefore,  for  this  gap  in  the  contour  caused 
by  the  zeroth  fringe,  the  effects  of  electrons  just  cancel  those  of  ions 
and  atoms.  A  calculation  may  be  made,  based  on  certain  reasonable 
assumptions,  which  will  set  an  upper  bound  on  the  electron  concentration 
at  this  point.  We  consider  only  the  l/l6-atm  case. 

Combining  the  results  for  refractivities  of  argon  and  of  free 

2 

electrons  given  by  Alpher  and  White  ,  and  assuming  that  ions  and  atoms 
have  equal  refractivities,  we  can  write 

pp  ji 

t^(p)  -  pCPq)  =  -  1.53  X  (10)  +  2.84  X  10  ^  (p  -  p^)/PgYP  (1) 

where  p  represents  refractive  index,  is  electron  concentration, 

o  and  p  are  densities  in  the  disturbed  and  undisturbed  beams 
o 

respectively,  p„„p  is  the  density  of  argon  at  standard  temperature  and 
bir  _  c 

pressure,  and  wavelength  X,  =  5*46  x  10  cm  has  been  used  in  evaluation 
of  the  constants. 

For  an  axisymraetric  shock  it  is  well  known  that  fringe  shift  is 
zero  at  t  ii0  sliOClC  but  rises  abruptly,  with  a  vertical  tangent  just 
Inside,  because  the  jiath  in  the  disturbance  increases  rapidly  while 
index  jump  remains  nearly  constant.  A  zero  shift  through  the  shock 
can  only  he  achieved  if  the  difference  in  in  indices  vanishes;  thus,  the 
left  side  of  Equation  (l)  is  zero  at  the  point  in  question. 
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The  head  shock,  SI,  is  seen  from  Table  I  to  have  M  si  5  there¬ 
fore  may  be  considered  a  strong  shock  wave.  The  density  jump  across 

such  a  shock  is  given  by  p^/p-,  =  (7  +  l)/(7  -  !)•  With  7  =  5/3  for 

^  X  jrj 

argon  and  with  p^  =  =  Pgqip/l^  immediately  find  =  4  x  10 

electrons/cc  for  the  1/16  atm  case.  This  electron  concentration  is  about 
the  same  as  that  of  the  first  contour  exhibited  in  W-  Further  back 
along  the  zeroth  fringe,  it  is  plausible  that  the  mean  electron  concen¬ 
tration  must  increase  continuously  to  counteract  the  positive  shift 

caused  by  the  positive  index  integrated  along  the  increasing  transverse 

17 

path  within  the  shock.  Thus  the  =  4  x  10  contour  shoiild  be  further 
and  further  above  the  zeroth  fringe.* 

Nearer  the  axis  and  just  beyond  the  expanding  wire  material,  the 

zeroth  fringe  appears  again,  this  time  as  a  narrow  band  bordered  on  the 

outside  by  the  6  =  -0.5  contour.  Here  the  interpretation  in  terms  of 

electron  densities  is  difficult  because  the  integrated  effects  of  the 

compressed  ions  and  neutrals  behind  the  shock  SI  together  with  electron 

effects  indicated  by  the  various  negative  fringe  contours  must  be 

included  to  yield  the  zeroth  fringe.  Lacking  a  quantitative  analysis, 

one  can  argue  that  the  abrupt  positive  jump  of  l/2  fringe  at  the  zeroth 

fringe  boundary  represents  a  compressional  shock  discontinuity.  This 

is  corroborated  by  the  sharply  hooked  fringes  of  the  corresponding 

multiple  fringe  interferograms  Pi-  The  appearance  of  compressional 

shock  SII,  argues  that  negative  shifts  caused  by  electrons  at  this 

boundary  are  negligible.  This  in  turn  implies  electron  densities 

W  ^  4  X  10^'^/cc. 
e  ' 

The  conclusion  could  be  false  if  SII  represents  a  density  jump 
Pg/Pl  >  4,  which  could  only  be  accounted  for  by  assuming  an  appreciable 
concentration  of  some  massive  atomic  or  molecular  species  behind  the 
shock,  e.g.  copper  atoms  or  argon  molecule  ions  Ag  The  density 


*This  is  easily  proven  if  shock,  electron  cloud  boxmdary  and  zeroth 
contoiar  occur  in  this  order  and  trace  similar  power-law  curves. 
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jump  is  probably  not  much  larger  than  ^4-  since  at  late  times  t  >  8  4isec, 
the  SII  boundary  is  broadened  and  nearly  erased  by  what  appears  to  be 
an  increase  in  electron  concentration.  This  effect  may  conceivably  be 
ascribed  to  increased  ionization  caused  by  the  relatively  weak  SII 
shock  as  it  passes  over  high  temperature  argon  atoms  already  in  excited 
states  produced  by  the  peripheral  arc.  The  multiple  fringe  inter- 
ferograms  especially  well  show  this  broadening  of  the  boundary  and 
the  accompanying  luminosity. 


5.2  Temperature  Behind  the  Head  Shock  wave. 


0.0  O.I-I  0.jr  i  It.  JH.  LH-iH-t- 1  IV  iipp  j 


as  may  be  further  corroborated  by  inspection  of  the  figures  in  [^1^  ; 
furthermore ,  the  argument  of  the  previous  paragraph  implies  that  SI  lies 


i  •nr*r*c»n  G  *1  n  crl '\r  on+‘ <=  r'l^rrh/^n-r*  rji -J- h  M  ~ 
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rsn+.oiiT’  with  M  =  4  X  10  /cc .  Since  the  shock 
e  ' 


is  compressional  and  strong  the  density  immediately  behind  it 

is  almost  k  p  .  Thus  for  the  l/l6  atm  case  the  argon  atom  density  is 
^  l3 

approximately  6  x  10  /cc  and  ionization  remaining  behind  the  shock 
should  be  less  than  7^*  The  shock  itself  is  too  weak  to  cause  further 
ionization  unless  appreciable  numbers  of  excited  argon  atoms  exist  in 
the  undisturbed  region.  We  discount  this  possibility  because  of  the 
absence  of  luminosity.  We  may  therefore  regard  SI  as  a  strong  shock 
wave  propagating  into  a  perfect  gas  of  7  =  5/5  at  room  temperature. 


With  these  assumptions,  the  temperature  ratio  across  a  strong 
shock  wave  traveling  at  Mach  M^^  is  given  by  T^/T^^  =  ^-y{y-l)/(y+l)^^ 

Substitution  shows  T^/T^^  slightly  less  than  8  and  T^  2400  °K. 

The  shock  wave  SI  is  approximately  straight  in  the  l/52  and  I/16 
atm  cases  and  therefore  proceeds  at  nearly  constant  Mach  number  over 
most  of  the  field  of  view.  Two  inferences  may  be  made:  l)  the  shock 
is  receiving  substantial  amounts  of  energy  from  the  expanding  interior 
regions;  otherwise,  it  would  decelerate  and  show  appreciable  curvature, 
and  21  the  firas  nassed  over  bv  +.Vie  shnek  .  ■miio-Vilv  'hnnnHe.rl  t  non  Ho  ■h^r  -l-Vio 

zeroth  fringe,  is  heated  at  least  to  the  temperature  Tg  ~  2400°K. 
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Estimate  from  Blast  Wave  Theory. 

Some  doubt  may  remain  whether  the  head  shock  wave  could  conceivably 
have  originated  at  the  expanding  wire  like  a  blast  wave  and  remained 
obscured  by  the  arc  plasma  until  appearing  beyond  the  zeroth  fringe  contour. 
The  measured  properties  of  the  wave  are  incompatible  with  such  a  view  as 
will  appear  from  the  following  argument. 

5 

From  the  equations  of  the  blast  wave  theory  for  cylindrical  shocks 
and  the  perfect  gas  law,  an  expression  for  temperature  ratio  across  the 
shock  can  be  found.  It  is 

m  .  V  . 

^  ,  Mnil  .  Mtil  (^Y  .4-  ,  (2) 

1  (r+i)  (r+i)  Pq  j  t 

where  the  symbols  have  their  usual  meaning,  t  is  time  and  the  velocity 

of  sound  in  the  undisturbed  medium.  The  constant  e  =  2.25  is  the  result 

4 

of  some  improved  calculations  that  furnish  values  appropriate  for  argon 
and  correct  certain  numerical  errors  in  Lin’s  tabulated  functions. 

The  first  equality  of  (2)  occurs  because  the  theory  assumes  the 
strong  shock  conditions  at  the  boundary,  and  therefore  must  agree  with 
the  result  of  if  3.2.  The  second  equality  expresses  the  ratio  in  terms 
of  the  time  and  the  constants  characterizing  the  assumed  shock  wave. 

r-T  ^  :z 

Using  E  =  0.5  X  10  ergs/cm  from  Table  II,  p  =  .105  x  10  gms/cc  and 

4  /  °  2 

a  =  3*2  x  10  cm/sec  the  quantity  corresponding  to  M-,  becomes 

o  2 

7l/t(psec).  Agreement  with  the  previous  value,  25,  is  achieved 

either  by  taking  t  'V'  3  psec  or  multiplying  E  by  factor  I/9.  Either 

alternative  contradicts  the  observations  and  must  be  considered 

unacceptable  as  a  true  description  of  the  shock  wave.  The  first  implies 

an  origin  in  time  2  psec  ahead  of  that  actually  observed.  The  second 

implies  a  shock  radius  R  only  l/5  of  that  observed,  through  the  dependence 
1  1 
If  2 

R  a  E  t  . 
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This  estimate  given  hy  blast  wave  theory  indicates  convincingly 
that  the  head  shock  SI  cannot  have  originated  on  the  axis  at  the  true, 
experimental  origin  in  time,  a  fact  we  have  already  appreciated  from 
a  different  point  of  view. 

$.4  Thermalizatlon  Times. 

It  is  of  great  interest  to  know  whether  conditions  inside  the 

peripheral  arc  approach  local  thermodynamic  equilibrium,  that  is  whether 

electrons,  ions  and  neutrals  are  closely  in  local  temperature  equilibrium 

or  not.  For,  if  so,  we  may  conclude  that  the  head  shock  wave  could  be 

produced  by  thermal  expansion  processes  of  the  usual  type.  If  not,  we 

may  then  consider  whether  the  shock  is  of  the  "electron  driven"  type 

5 

reported  by  Fowler  et.al  where  the  hot,  expanding  electron  gas  draws 
the  cold  ions  and  neutrals  with  it  by  means  of  Coulomb  and  collisional 
forces. 

In  this  section  we  employ  known  experimental  data  and  theoretical 
expressions  to  calculate  critical  time  intervals  for  approach  to  .thermal 
equilibrium  between  electrons,  neutrals  and  ions  for  various  assumed 
collision  processes. 

We  consider  first  the  conditions  which  may  exist  at  the  edges  of 
the  peripheral  arc  and  in  the  region  where  the  head  shock,  SI,  appears 
to  originate.  Our  experimental  data  M  and  previous  considerations 
indicate  that  the  electron  density  is  sufficiently  low  so  that  electron- 
atom  collisions  may  be  considered  to  predominate.  Even  though  electron- 
electron  collisions  may  occur  only  once  for  every  15  electron-atom 
collisions,  because  of  the  great  disparity  of  masses  (2m/M  *'-'l/40,0CX)) , 
the  electrons  lose  little  energy  to  the  atoms,  yet  remain  in  temperature 
equilibrium  themselves. 

Making  lise  of  the  definitions  and  data  presented  by  Brown^,  we 
present  in  Table  III  the  results  of  calculations  for  electron-atom 
collisions  in  argon.  Here  V  is  electron  energy  in  ev,  T^  is  electron 
temperature  in  °K,  v  is  electron  velocity  in  cm/sec,  the  average 
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number  of  electron-atom  collisions  per  cm  of  electron  travel  per  mm 

Hg  pressure  at  0°C,  v  is  ih®  electron-atom  collision  frequency, 

p  =  275  p/T  where  p  =  47*5  mm  (1/16  atm),  X  =  v/v  is  the  electron 
°  4  /  e  c 

mean- free -path  and  r  =  4  x  10  /v^  is  the  average  time  during  which  an 

electron  would  transfer  its  energy  to  atoms  i,e.  it  is  assumed  to  take 

in  the  neighborhood  of  40,000  electron-atom  collisions  to  bring  electrons 

into  thermal  equilibrium  with  atoms. 

TABLE  III 
1/16  ATM  ARGON 

1 


v2 

T 

e 

(°K) 

V 

(cm/sec) 

P 

c 

V 

c 

coll/ sec 

X 

e 

(cm) 

T 

(psec) 

1 

8  X  10^ 

Q 

0.6  X  10 

4 

0.1  X  10^^ 

6  X  10"^ 

4.0 

2 

32  "  " 

1  ^2  "  " 

27 

1.4  "  " 

9  X  lO"^ 

0.3 

5 

72  "  " 

1.8  "  " 

72 

5.6  "  " 

J  tt  11 

0.07 

5.5 

96  "  " 

2..1  "  " 

83 

8.0  "  " 

2.6"  " 

0.05 

Because  of  the  Rarasauer  effect  in  argon,  the  probability  of 
collision  varies  by  a  factor  of  more  than  20  as  electron  energy 
goes  from  1  to  12.  Inspection  of  the  last  column  indicates  a  hundred¬ 
fold  variation  in  time  elapsed  to  obtain  thermal  equilibrium  between 
electrons  and  atoms ,  One  may  not  conclude  that  the  longer  thermalization 
times  prevail,  because  the  axial  electric  field  may  be  of  the  order  of 

several  kv/cm  during  the  formation  of  the  peripheral  arc  and  for  several 

7 

hundred  nanoseconds  thereafter  .  Since  it  seems  unlikely  that  the 
average  electron  energy  is  less  than  2-3  ev,  we  may  tentatively  conclude 
that  the  thermalization  time  is  certainly  not  longer  than  the  time  for 
appearance  of  the  head  shock  wave  and  that  the  head  shock  is  probably 
mainly  thermal  in  origin. 

As  for  conditions  near  the  center  of  the  peripheral  arc  and  in 
regions  where  the  ionization  exceeds  10^,  we  have  recourse  to  theoretical 
expressions  derived  for  completely  ionized  plasmas  in  which  collision 
processes  are  conditioned  on  long-range  electrical  forces. 
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g 

Making  use  of  Spitzer’s  Equations  (5  -  26)  and  (5  -  51)  we 

calculate  the  electron,  self-collision  time  t^  and  the  equipartition 

time  t  ,  between  electrons  and  argon  ions.  Because  of  the  massiveness 
eq 

of  argon  ions  relative  to  electrons  and  less  importantly  because  of  the 
relatively  low  temperature  we  assume  for  the  ions  initially  (T-,  'v-  500°K) , 


each  expression,  t^  or  t^g^,  is  proportional  to 


5/2 


and  thus  increases 

rapidly  with  electron  temperature.  In  order  to  establish  what  may  be 

an  upper  bound  for  these  times,  we  assume  the  mean  electron  energy  to 

5  o 

be  about  10  ev  (t^ 10"^  K).  Then  taking  the  number  density  of  electrons 

l8  12  8 

or  ions  to  be  10  /cc,  we  obtain  t  "^10  sec  and  t  5  x  10  sec. 

c  eq 

The  self -collision  time  for  argon  ions  at  room  temperature  is  larger  than 


2  x/p 

that  for  electrons  in  the  ratio  (40  x  1823)  (3  x  10  )  '  and  as  the 

temperature  of  the  ions  rises  toward  10^  °K,  the  second  factor  would 
approach  unity.  Thus  at  the  outset  the  self -collision  time  for  argon 

ions  is  actually  .04  that  of  electrons  but  rises  to  27O  time  as  much. 

-9 

It  never  exceeds  10  sec. 


These  calculations  support  the  view  that  initially  the  electrons 
and  ions  are  in  temperature  equilibrium  among  themselves  and  come  into 
temperature  equilibrium  with  each  other  in  not  longer  than  50  nanosec, 
a  time  represented  by  about  1/5  of  the  first  light  fringe  at  the  tip 
of  the  l/16-atm  interf erogram,  Figure  2.  Thus  we  shoiild  expect  conditions 
within  the  peripheral  arc  closely  to  approximate  local  thermodynamic 
equilibrium  until  the  electron  density  falls  to  a  value,  arbitrarily  taken 
here  as  10^  ionization,  below  which  plasma  behaviour  can  no  longer  be 
expected.  Even  in  the  outer  regions  of  low  ionization,  thermalization 
times  for  the  more  energetic  electrons  are  of  the  order  of  50  -  100 
nanosec,  as  may  be  seen  from  Table  III. 

3.5  One -Dimensional  Unsteady  Flow. 

The  fact  that  local  thermodynamic  equilibrium  can  be  expected  from 
point  to  point  throughout  most  of  the  peripheral  arc  region  suggests 
another  approach  to  estimating  the  average  temperature  there. 
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Let  us  assume  that  just  prior  to  the  emergence  of  the  shock  SI, 
the  annulus  of  the  peripheral  arc,  which  may  be  considered  as  nearly 
a  right  circular  cylinder,  contains  argon  gas  at  high  temperature  and 
pressure  but  at  or  near  the  original  undisturbed  density.  This 
situation  could  be  realized  because  of  the  rapidity  of  heating  and  the 
massiveness  of  the  gas.  Then,  in  the  next  interval  of  time,  a  few 
tenths  psec,  expansion  of  this  heated  gas  gives  rise  to  the  head  shock 
wave.  If  the  outer  radius  is  sufficiently  large  and  the  time  interval 
short  enough  so  that  the  expansion  wave  propagates  only  a  small  distance 
back  into  the  high  temperature  reservoir,  the  results  of  unsteady  flow 
theory  applicable  to  the  one-dimensional  shock  tube  may  be  used. 

9 

Resler  et.al.  developed  an  expression  for  the  pressure  ratio 
Pj^Pf  between  reser’voir  pressure  Pj^  and  pressure  in  the  'undisturbed 
region  ahead  of  the  shock.  Their  Equation  (7)  is  a  function  of  the 
shock  Mach  number  the  velocities  of  sound  a^,  a|^  and  the  gas  constants 
^1’  Under  the  assumption  already  m^de  that  and  =  5/5 

with  the  perfect  gas  law  and  a  =  (7  RT)  ,  we  can  obtain  an  implicit 
expression  for  the  temperature  ratio  a  function  of  M^. 

A  numerical  solution  of  this  equation  shows  that  'v  5  corresponds  to 
T|^/t^  Thus,  for  our  l/l6-atm  case  the  temperature  of  the  heated 

plasma,  before  appreciable  expansion,  is  indicated  to  be  about  22,000  °K. 
It  is  interesting  to  observe  that  this  is  the  temperature  for  which 
Olsen^'^  finds  a  maximium  in  the  conductivity  of  a  thermalized  argon  plasma 
at  one  atmosphere  pressure.  He  observes  temperatures  of  this  order  near 
the  anode.  Increasing  total  current  has  the  effect  of  enlarging  the  arc 
and  the  high  temperature  (20,000  °K)  region  near  the  anode  but  does  not 
increase  the  maximum  temperature  appreciably. 

5-6  Sound  Speed  in  the  Arc  Region. 

Close  examination  of  the  interferograms  for  the  I/16-  and  l/52-atm 
cases  shows,  that  the  weak  wave  WS  originates  on  the  expanding 
wire  surface  at  the  time  that  the  inward  expanding  arc  boundary  intersects 
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the  wire.  Thence  it  expands  along  a  nearly  straight  path  which  crosses 
the  S  =  -  0.5  fringe  contour  and  shows  a  pronounced  shift  toward  larger 
radii  on  the  downstream,  or  later  time,  side.  Comparison  shows  that 
the  fringe  shift  across  WS  is  negative  and  suggests  that  a  detectable 
Increase  in  electron  density  has  occurred.  The  disturbance  is  not 
easily  traced  into  the  compressional  region  behind  SI  probably  because 
the  compressional  density  change  is  too  weak  to  detect. 

Comparison  between  the  l/l6-atffl  and  higher  pressure  interferograms 
shows  that  the  flare  in  the  first  negative  fringe,  with  attendant 
negative  fringe  shift,  can  be  discerned  in  several  cases.  This 
observation  suggests  that  similar  disturbances  are  present  in  these 
also,  but  cannot  be  seen  at  early  times  because  at  higher  pressures 
the  arc  region  is  much  smaller  and  because  luminosity  obscures  the 
field  of  view. 

The  cause  of  disturbance  WS  is  not  known  and  not  readily  inferred. 
Because  it  originates  from  the  most  rapid  arc  process  and  at  early 
enough  times  so  that  complete  thermalization  may  not  have  occurred, 
the  suspicion  is  aroused  that  it  may  be  a  shock  or  abrupt  disturbance 
of  the  "electron  driven"  category  described  by  Fowler  et.al^. 

The  facts  that  the  disturbance  WS  has  a  moderate  velocity,  reflects 
from  the  dense,  metal  boundary  and  causes  a  slight  but  detectable 
increase  in  electron  density  suggests  that  it  is  hydrodynaraical  in 
nature;  for  a  weak  shock  in  a  moving,  high-temperature  gas  or  plasma 
could  propagate  at  several  times  the  speed  of  sound  in  the  eunbient  cold 
gas  and  by  compression  could  slightly  increase  the  temperature  and  the 
accompanying  ionization  in  the  gas  over  which  it  passes. 

We  shall  assxime  that  WS  is  indeed  a  weak  distiirbance  and  that  it 
propagates  in  a  high  temperature  plasma  which  is  in  local  thermodynamic 
equilibriiun.  Its  velocity  relative  to  the  fluid  in  which  it  moves  is 
then  slightly  greater  than,  or  equal  to,  the  wave  velocity  in  the  medium. 
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According  to  the  treatment  of  positive-ion  oscillations  which  preserve 
electrical  neutrality,  given  by  Spitzer^^,  this  velocity,  for  the  case 
of  argon  plasma  is  given  by  a^  =  (2y'kn!^) /m^  where  m^  =  m^  is  the  tonic 
or  atomic  mass. 


The  sound  speed  in  the  undisturbed  region  is  just  a  =  (yhT  )/ni.  i 

p  O  O  A 

thus  we  find  2T  /T  =  (a  /a  )  .  Taking  a  to  be  the  velocity  of  the 
p'  o  '  p'  o  P 

weak  shock,  WS,  from  Table  I  and  a^  =  0.52  ram/psec,  we  obtain 

T^/t^  =  86  or  T^  26,000  K.  The  agreement  with  our  previous  estimate, 

obtained  from  one-dimensional  unsteady  flow  theory,  seems  quite 

satisfactory. 


We  have  not  made  any  allowance  for  the  fact  that  the  fluid  supporting 
WS  may  be  in  motion.  We  may  infer  that  it  is  in  a  state  or  practically 
uniform  motion  from  the  straightness  of  WS  and  the  absence  of  other 
discontinuities . 

The  inner  boundary,  SII,  is  moving  at  1  mm/psec  and  the  outer 
boundary,  SI,  at  1.5  mm/usec.  If  SI  had  passed  over  the  entire  region 
traversed  by  WS  the  velocity  of  the  fluid,  given  by  the  strong  shock 
relations  would  be  2/ (7  +  l)  =  3A  times  the  velocity  of  the  shock. 


It  seems  more  likely  that  only  the  latter  end  of  WS  could  be 
affected  by  SI  and  that  the  fluid  traversed  by  WS  is  practically  at 
rest.  Nevertheless,  the  above  argument  makes  it  appear  quite  certain 
that  the  fluid  velocity  in  the  region  of  WS  cannot  exceed  1.1  mm/psec. 

If  we  make  this  correction  to  the  velocity  of  WS  the  temperat\ire 

estimate  becomes  T  /T  =46  and  T  ‘^14,000  °K.  This  value  may  be 
p  o  P  ' 

regarded  as  somewhat  too  low. 


The  possibility  that  WS  is  a  magnetohydrodynamic  disturbance, 
dependent  on  magnetic  field,  cannot  be  ruled  out  but  seems  unlikely 
for  the  following  reasons.  WS  is  practically  straight  throu^out  a 
time  interval  long  enough  for  at  least  one  current  (and  field)  reversal 
to  take  place.  No  corresponding  variation  in  slope  of  WS  can  be  seen. 
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Since  the  current  has  transferred  from  the  wire  to  the  arc  and  the 


current  density  may  he  considered  reasonably  uniform,  the  magnetic 
field  will  rise  linearly  with  radius  through  much  of  the  region  traversed 
by  WS.  Strong  dependence  on  magnetic  field  would  argue  changing  slope 
of  WS  which  is  not  observed.  We  conclude  that  WS  can  be  at  most  only 
weakly  dependent  on  magnetic  field  and  are  therefore  justified  in 
considering  it  to  be  an  acoustic  wave  in  the  plasma. 


^.7  Electron  Driven  Shock  Waves. 

According  to  Fowler  et.al.^  the  velocity  of  the  fast  luminous 
front,  which  they  indent if y  as  electron  driven,  is  given  by  V  ~  k  T^/M, 
where  Is  the  electron  temperature  and  M  is  the  ion  mass.  A  calcu¬ 
lation  with  Tg  '^2.5  ev,  which  corresponds  to  our  previous  temperature 
estimates,  gives  a  velocity  of  2  mm/psec.  This  is  comparable  with  the 
velocity  of  WS  but  considerably  smaller  than  the  lower  bound  on  average 
velocity  of  P  given  in  Table  I.  This  fact  suggests  that  the  effective 
electron  temperature  just  inside  the  front  surface  P  may  be  considerably 
larger  than  the  value  of  2.5  ev  which  ought  to  hold  in  the  arc  region 
where  ions  and  electrons  are  in  temperature  equilibrium. 

We  have  no  independent  measure  of  electron  temperature  just  behind 
the  arc  front  P,  but  setting  V  =  2  cm/psec  as  a  lower  bound  yields  an 
electron  temperature  T^  '^200  ev  in  the  front.  Such  a  high  value 
forces  us  to  reconsider  the  estimated  collision  and  equipartition  times 
calculated  earlier. 

The  interferograms  and  electron  contours  of  Cl]  show  the  leading 
edge  of  the  peripheral  arc,  in  which  ionization  rises  from  zero  to  10^, 
to  be  somewhat  less  than  l/5  mm  in  linear  dimension  and  when  extrapolated 
to  the  tip  about  100  nsec  in  duration.  In  this  narrow  border,  electron- 
atom  collisions  will  predominate  and  calculations  like  those  of  Table  III 
may  be  appropriate. 

If  we  extend  the  collision  cross-section  curve  for  argon^  with  an 
hyperbolic  approximation  which  fits  the  high-energy  tail  quite  well,  we 
find  that  the  time  t  for  electrons  and  atoms  to  come  to  equilibrium 
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approaches  a  constant  value  of  about  75  nsec  but  remains  in  the  range 
50  -  80  nsec  for  electron  energies  down  to  9  ev.  These  estimates  show 
clearly  that  as  ionization  rises  through  10^  in  the  arc  front,  all 
electrons  with  energies  higher  than  9  ev  can  approach  temperature 
equilibrium  with  atoms.  Those  initially  at  200  ev  will  fall  rapidly 
in  temperature  in  the  same  interval  unless  energy  is  continuously  added. 
Because  the  relative  number  of  electrons  is  small  the  resultant  mixing 
temperature  will  be  much  lower  than  the  maximum  electron  temperature. 
Beyond  the  10^  ionization  contour,  plasma  conditions  will  begin  to  . 
prevail  and  the  estimates  of  <r  may  be  used.  For  T^  10  ev  the 

equipartition  time  is  less  than  50  nsec;  therefore,  in  a  total  time  not 
exceeding  I50  nsec  we  would  expect  the  high- temperature  electrons  to  have 
cooled  and  come  to  thermal  equilibrium  in  the  arc  plasma.  Our  previous 
estimates  put  the  equilibrium  temperature  at  about  2.5  ev. 

To  summarize,  if  we  regard  WS  as  the  reflection  of  P  from  the  wire 
boundary,  then  P  is  in  some  ways  like  a  fluid  mechanical  shock  wave. 
Fiorthermore ,  the  velocity  of  P  and  knowledge  of  electron  densities  behind 
P  indicate  that  in  the  leading  edge  of  the  arc,  less  than  0.1  psec  wide 
electron  temperatures  are  of  order  200  ev  and  are  far  higher  than  those 
for  the  surrounding  atoms  and  ions  which  were  initially  at  room 
temperature.  We  conclude  on  the  basis  of  the  foregoing  argument  that 
the  assumption  that  the  arc  front  P  is  an  electron  driven  shock  wave  of 
the  type  discussed  by  Fowler  and  co-authors  is  compatible  with  the 
experimental  data  and  offers  a  consistent  explanation  of  the  leading 
edge  of  the  arc. 
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4.  SUMMAJ^Y 


As  a  sequal  to  our  earlier  interferometric  study  we  present  single 
fringe  interfere grams  of  exploding  copper  wires  at  ambient  pressures  of 
1/8,  l/l6,  1/32,  atm  of  argon.  Features  of  the  disturbance  not  easily 
recognized  in  multiple  fringe  pictures  stand  out  clearly.  In  particular 
it  is  possible  to  distinguish  the  peripheral  arc  boundary  P,  the 
compressive  head  shock  wave  SI,  a  fast  weak  disturbance  WS  and  a 
compressive  shock  at  the  expanding  boundary  of  the  copper  metal  SII. 
Representative  velocity  values  obtained  from  the  l/l6-atm  case  are 
used  in  a  semi-quantitative  analysis. 

The  following  main  points  are  noted: 

a.  The  fluid  disturbance  can  be  divided  roughly  into  three 
regions;  a  compression  bounded  on  the  outside  by  the  head  shock  and 
internally  by  the  zeroth  fringe  contour;  an  annular  plasma  bounded  by 
the  zeroth  fringe  contour  and  internally  by  the  expanding  wire;  and 

the  expanding  wire  vapor,  itself  the  nearest  of  the  regions  to  the  axis. 

b.  Strong- shock-wave,  boundary  conditions  indicate  that 
temperature  behind  the  head  shock  is  about  2500  °K.  Application  of 
blast  wave  theory  confirms  that  the  flow  is  still  gaining  energy  during 
the  time  of  observation  and  that  it  originates  not  at  the  wire  surface 
but  on  an  annulus  about  1  cm  in  diameter  formed  by  the  expanding 
peripheral  arc.  Ionization  rises  from  a  few  percent  just  outside  the 
zeroth  fringe  contour  to  values  in  the  range  50  -  70^  near  the  center 
of  the  arc, 

c.  Estimates  of  thermalization  times  for  electron-atom  and 
electron- ion  collisions  show  that  local  thermal  equilibrixmi  should  be 
well  approximated  everywhere  except  near  the  arc  boundary  P. 

d.  Assumption  that  the  head  shock  is  Initiated  by  a  one¬ 
dimensional  adiabatic  expansion  of  a  hot  cylindrical  plasma  formed  by 
the  peripheral  arc,  leads  to  an  estimated  mean  temperature  of  the  plasma 


22 


of  about  25,000  °K  (2.5  ev).  Identification  of  the  weak  shock  wave  WS 
with  a  sound  speed  disturbance  of  the  plasma  yields  almost  the  same 
mean  temperature  in  the  plasma. 

e.  Finally,  consideration  of  conditions  immediately  behind 
the  leading  arc  boundary  P,  indicates  that  thermal  equilibrium  between 
electrons  and  ions  does  not  hold,  that  electron  temperatures  as  high 
as  200  ev  may  occur,  and  that  the  arc  boundary  may  reasonably  be  thought 
of  as  a  shock  wave  of  the  electron  driven  type. 

F.  D.  BENNETT 
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Fig.  1  alngle  Fringe  StreaJc  Interferograra. 

4  rail  Cu  wire  at  20  kv,  l/52-atra  argon. 
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Single  Fringe  Streak  interferogram. 

4-mil  Cu.  Wire  at  20  kv,  l/l6-atm  argon. 


TTt  cf  O 

I- 


Fig.  3 


Single  Fringe  Streak  Interf erogram. 
k-mil  Cu  wire  at  20  kv,  l/8-atni  argon. 
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